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CLAIMS 

An electromagnetic simulation algorithm, for 
determining the electromagnetic wave scattered by 
a body in a monof requency situation, from a 
meshing of said body and from the electromagnetic 
excitation, 

characterized in that it comprises at least: 

(a) a determination of a matrix M, the so-called 
interaction matrix, whose coefficients are 
determined from the meshing of said body; 

(b) a determination of a precondi tioner Z of the 
matrix M, this precondi tioner being the 
matrix formulation of the operator 'J*9A.J 
where ftyt is the operator associated with M, 
J the operator vector product with the 
normal to the surface of said body and J* the 
operator adjoined to J ; 

(b) a determination of the currents which flow 
around the surface of said body, through an 
iterative algorithm of conjugate gradient 
type using said precondi tioner Z, the 
iterative algorithm making it possible to 
solve an equation, the so-called boundary 
integral equation of electromagnetism, 
written in matrix form in the following 
manner : 

MU = L 

where U is a vector which we seek to 
determine, whose coefficients represent the 
surface currents, and L is a known vector, 
whose coefficients represent the 

electromagnetic excitation ; 

(c) a determination of the wave scattered by said 
body, from said surface currents. 

The electromagnetic simulation algorithm as 
claimed in claim 1, characterized in that the 
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coefficients of the vector U are expressed in the 
usual basis of the Ravi art -Thomas space. 

The electromagnetic simulation algorithm as 
claimed in one of the preceding claims, 
characterized in that the preconditioner Z is 
determined implicitly 

The electromagnetic simulation algorithm as 
claimed in one of the preceding claims, 
characterized in that the preconditioner Z is 
defined by the following relation: 
Z= C JMJ 

where J is a matrix formulation of the operator J 
and t J the transposed matrix of J. 

The electromagnetic simulation algorithm as 
claimed in claim 4, characterized in that the 
matrix J is defined by the following relation: 

J = (1 - M6 {o Ms] U M4JU2 oj loj 

where 
M1,=Jp,(A)-<p,(A)dS A 
M2, = Jv ( / J (A)div A (^)dS A 
M3, = jfotA) a z(A))- <p,(A) dS A 
M4,=Jv j (A)9,(A)dS A 

M5,= Je,(A)e 1 (A)ds A 

where 



(Pi is the usual basis of the Raviart-Thomas space; 



Gi is the basis of the trianglewise affine 
functions ; 



is the basis of the trianglewise constant 
functions ; 

A a point belonging to the surface T of said body; 

z(A) is a vector with unit norm, normal to the 
surface of said body at the point A, and oriented 
outward . 

The electromagnetic simulation algorithm as 
claimed in one of the preceding claims, 
characterized in that the iterative algorithm used 
is a fast algorithm, of the multilevel multipole 
method type. 

The electromagnetic simulation algorithm according 
to one of claims 1 to 5 , characterized in that the 
iterative algorithm used is a fast algorithm, of 
the adaptive integral method method type. 

The electromagnetic simulation algorithm as 
claimed in one of the preceding claims, 
characterized in that the body is an antenna for 
which one seeks to determine an optimal shape, by 
using the simulation algorithm in an antenna 
design tool. 



The electromagnetic simulation algorithm as 
claimed in one of claims 1 to 7 , characterized in 
that the body is an object of known shape for 
which one seeks to determine the radar cross 
section (RCS) . 



